

/ 

(■> 

’1 

ll 

rH 

fi in 

VO 

=- * 

BbhhTD 

11! II 

H 

• 

NATIONAL ADVISORY COMMITTEE 

~ k,- = 

i 

FOR AERONAUTICS 

V 


TECHNICAL NOTE 


co 
M- 
ro 
— J 


NO. 1651 

SUPERSONIC NOZZLE DESIGN 
By J. Conrad Crown 


I 

r* 

k 


Ames Aeronautical laboratory 
Moffett Field, Calif. 




WASHINGTON 

joke 1948 



U l-Vn 1 l! b . > ’ 


n q ^ n 

, ' 2 jJj iJti 11 


'• (T- - , <1 

4S1L1 lU U L 


TECH LIBRARY KAFB, NM 


ifcCH LibmARY KAt-fcJ, NWl 


NATIONAL .ADVISORY COMMITTEE EOS AERONAUTICS 


MM 1 


TECHNICAL NOTE NO . ^ I65I 


SUPERSONIC NOZZLE DESIGN 
By J. Conrad Crown 


SUMMARY 

The theory of supersonic flow in nozzles is discussed, emphasis 
heing planed on the physical rather than the mathematical point of 
view. Methods, "both graphic and analytic, for design in g nozzles are 
described together with a discussion of design factors. In addition, 
the analysis of given nozzle shapes to determine velocity distribution 
and possible existence of shock waves is considered. A description of 
a supersonic protractor is included in conjunction with a discussion 
of its application to nozzle analysis and design. 


INTRODUCTION 

One of the major problems in the design of a supersonic wind 
tunnel is the determination of the contours of the supersonic nozzle 
so that parallel and uni form flow in the test section may be assured. 
Consequently, it is not surprising that the literature contains 
numerous papers on the subject of supersonic nozzle design. These 
vary widely in their degree of complexity and general availability. 

It Is the purpose of this report to discuss these various methods 
and present a guide for nozzle design. Only two-dimensional nozzles 
will be considered. 

The most prominent method for determining nozzle contours is, 
perhaps, that of Prandtl and Busenann (reference l). The usual 
presentation of their method of characteristics is rather mathematical 
in nature. (See,e.g., Preiswerk, reference 2 .) In order to provide 
the designer with a clearer physical picture of the flow in a nozzle, 
a different interpretation of the Prandtl— Bus emann method is presented. 
The diverse systems for constructing nozzle shapes by this method are 
also presented, together with certain ramifications and supplementary 
useful information. 

The Eoelsch method (reference 3) is included because its analytic 
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nature offers certain advantages.. : These will he discussed later. 
Shapiro (reference 4 ) has still another approach to the problem. His 
method, due to its approximate nature and because it -has no special 
advantages, will not be considered. 


BASIC THEORY 

It is well known that. In a purely contracting flow, the maximum 
unif orm velocity that can be achieved across any section is that 
corresponding to the local velocity of sound. Further increases in 
velocity can be obtained only by subsequent expansion of the Btream. 

The essential and relevant features of a channel designed to 
produce supersonic flow are shown in figure 1 . A compressible fluid 
at virtually zero velocity in the settling chamber is accelerated 
through the contraction section to sonic speed in the throat where, 
if the contraction section is properly designed, the flow is uniform 
and parallel. The fluid is then expanded in the nozzle until the 
desired Mach number is reached in the test section where the flow 
is again uniform and parallel. In the analysis, the nozzle itself 
is divided at the inflection point of the wall into two sections: 
initial and terminal. 

It should be noted that there is one additional prerequisite 
for the establishment and perpetuation of supersonic flow. This is 
the maintenance of at least the mi ni mum pressure ratio between the 
settling chamber (pressure = p Q ) and the test section (pressure = pt) 
from reference 5, page 2 6 

I- + 

where Mj. is the Mach number in the test 
ratio of the specific heats of the gas. 

An irrotational, nonviscous supersonic flow through a two- 
dimensional nozzle may be treated by means of a few simple considera- 
tions. First, consider an incident uni dimensional supersonic flow over 
a single curved surface. The change in local Mach number between 
any two points is a function only of the change in direction of the 
stream between the points or the change in direction of the tangents 
to the surface at these given points. ' To consider the flow field 
between two curved surfaces, however, it is convenient to replace 
each surface by an infinite number of infinitesimally long straight- 
line segments, or a finite number with discreet but small length. 


y 

y-i 

(1) 

section and 7 is the 
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Each adjacent pair of lines thus constituted forms a comer. The 
supersonic flow about a corner is a classical problem and its solution 
is known. The flow between two curved surfaces thus reduces to the 
determination of the combined effect of two sets of comers. This 
Introduces the problems of intersection and reflection of influence, 
or disturbance lines. In addition, the condition requiring uniform 
and parallel flow in the test section leads to the concept of 
neutralization of disturbance lines. The following sections will 
elucidate upon these concepts. 

Flow About a Corner 

The flow about a convex comer formed by two intersecting straight 
lines has been treated analytically by Prandtl and Meyer (reference 6, 
pp„ 2^3 -Qk6). For any such configuration, three regions of flow exist. 
These are indicated in figure 2. The flow is uniform and parallel up- 
stream and downstream of the corner in the regions I and HI bounded 
by the surface and the corresponding Mach lines. In the region II 
between these Mach lines, flow parameters are constant along radial 
lines (each of which is a Mach line) emanating from the vertex of the 
corner. 

The fundamental equation of flcrw-about— a-comer is (fig. 3) 

v . K taa-1 (s%2)- (90° - a) (2) 


where V is the expansion angl e or the angle through which the flow 
is turned in accelerating from a local Mach number of unity to any 
given Mach number M, a is the corresponding Mach angle 

. -l 1 
a = sin — 

M 


and 



Obviously, if V is known in any region, the Mach number is deter- 
mined by equation (2) and can be found. 

Let the subscripts 1 and 2 refer to conditions in regions 
I and III, respectively, of figure 2, then the angle through which 
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the flow is turned in accelerating from a Mach number Mi to M 2 , 
that is, in going from region I to HI, is 

8 = V 2 - V x (3) 

In other words, the change in expansion angle is equal to the absolute 
■value of the change in stream deflection through an expansion region 
due to a single comer. 

If the stream deflection angle 8 is small, then all the expan- 
sion may be considered to take place along the average Mach line as 
shown in figure This line, no longer a line of propagation of an 
infinitesimal disturbance, now takes tm certain characteristics of a 
shock wave; namely, the flow through it suffers a finite change in 
direction and Mach number. It is usually ref erred to as an expansion 
wave. Little error is introduced by making these assumptions and, as 
5 approaches zero, the error vanishes. It is convenient to define 
the strength of a wave as the angular deflection of the stream that it 
produces. This is numerically single valued for expansion waves and 
weak oblique shock waves. 


Flow Parameters 


Flow conditions are completely determined by the parameters v, 
the expansion angle, and 0, the stream angle relative to same datum 
line usually taken as the flow direction in the throat. These coordi- 
nates are usually written v , 0, or 



Intersection of Expansion Waves 

The problem of the interaction of the expansion waves from two 
opposed convex surfaces, such as the initial portion of a nozzle, 
may be considered in its elementary form: the intersection of two 

expansion waves as depicted in figure 5» 

It follows from reference 2, pp. 55—58* that the angular change 
in direction of the stream through an expansion wave is constant along 
its length regardless of the direction or velocity of the flow in front 
of the wave; that is to say that the expansion waves pass through each 
other mutually unaffected in strength, although their inclination is 
altered. Their effect on the flow may be determined by superposition 
of individual effects. 

Consider the two expansion waves shown in figure 5* S’ 0 ** con- 

venience, they are designated (l) and (2) and have strengths of +e 
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and —5, respectively. The upper streamline shown is deflected up 
through, an angle e by (l) and down through an angle & by (2) . The 
total angle through which it is deflected is thus +(e— 5) . Similarly, 
the lower streamline is deflected first downward by (2) then upward 
by (l) . Its final angle is the same as for the upper streamline and 
is equal to 6 + £- 5. In a like manner, the final expansion angle can 
be found to be increased by £+& for both streamlines . 


Reflection of Expansion Waves by a Wall 

CondLitions resulting from the reflection of an expansion wave by 
a bo undar y may be determined by utilizing the well-known mirror— image 
concept. Thus, the wall may be replaced by a stre aml ine in a fictitious 
flow comprised of the original flow, plus an image flow field, as shown 
in figure 6. The problem of the reflection of expansion waves by a 
wall then becomes that of the intersection of expansion waves. The 
latter problem was the subject of the preceding section. 

This concept may be applied in a converse manner in the design 
of symmetrical nozzles. In this case, the straight center line of the 
nozzle is replaced with a wall. Thus, the amount of work is halved. 


neutralization of Expansion Waves 

If a shock wave of infinitesimal strength is superimposed on an 
expansion wave of equal strength, (and by definition opposite in sign) , 
the flow is unchanged after passing through both. This is also very 
nearly true if the waves have a finite but small strength. Therefore, 
if at the point where an expansion wave hits the wall a compression 
wave of equal strength is created, the expansion wave will be 
neutralized. Such a compression wave can be created, as illustrated 
in figure 7j by an. angular change in direction of the wall equal to 
the strength of the given expansion wave. The direction of the 
deflection should be such as to form a concave comer. 


Flow in a Wozzle 

r T I hR flew throughout a two-dimensional nozzle can be determined 
by use of the previously discussed concepts. The flow coordinates 
in the nozzles shown in figures 8 and 9 aro presented to i l l u strate 
the method. While symmetrical nozzles are discussed predominantly 
herein, the concepts involved apply to supersonic flows in general. 

The angle between the wall at its inflection point and the 
center line (for symmetrical nozzles) is of importance in nozzle 
design. For shapes simulated by straight-line segments, this 
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inflection point appears as a region. Let the subscript 1 refer to 
conditions in this region immediately preceding the point at which 
neutralization first takes place. These positions are denoted by 
arrows in figures 8 and 9 * 

The following relation then becomes apparent: 


Vx + 01 = Vfc (4) 

In addition, it is obvious that the maximum value that 9 X can have 
for shock— free flow occurs when 6 1 = v x or 

^max = 2 " v t (5) 

It should be noted thAt, if the initial curve is not approximated 
by straight-line segments, 0 i can equal v x only for a nozzle which 
has an abrupt expansion at the throat as shown in figure 10. However, 
for such a nozzle, it is still possible for 9j_ to be less than V^-/ 2 , 
provided that some of the expansion waves are allowed to be reflected 
before they are neutralized. 

For any smooth initial curve, that is, with no discontinuity in 
ordinate or slope from the sonic section to the inflection point, 

V is greater than |0 | for 9^0. This condition appears to be 
violated in the nozzles shown in figures 8 and 9 > wherein there exist 
certain regions along the mil, where V equals 9. The explanation 
of this lies in the fact that the wall was simulated by a finite number 
of comers. The error introduced by this assumption is approximately 
given by 

0<V .-V <6 

exact approx 

where 5 is the angular deviation of each comer. In the cases 
illustrated in figures 8 and 9 j 5 equals 2°. Consequently, V 
is actually greater than 9 . This error is usually small and can be 
ignored without serious consequences. 

For any given Mach number, while there are an infinite number 
of satisfactory nozzles, there is one invariant parameter: the ratio 

of the areas of the test section and throat (reference 5 , 3^0 


A t 

A* 


1 [2 + (7-1) M t 2 -] 

Mfc L 


x_ 

2 



7+1 
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where A is cross-section area (or height in a two— dimensional, nozzle) , 
the * refers to conditions in the throat (sonic section), and the sub- 
script t refers to conditions in the test section. 


METHODS OF NOZZLE DESIGN AND ANALYSIS 
Busemann's Method 

Busemann’s method for designing nozzles (reference 2) consists of 
assuming an initial curve and finding the terminal curve required to 
give uniform and parallel flow in the test section at the desired Mach 
number . 

Tn order to design a nozzle for a Mach number M-fc, first find the 
corresponding expansion angle V-5. Assume an initial curve, and simu- 
late it by a series of preferably equiangular corners. Then, starting 
at the throat and proceeding downstream, determine the flow field in 
terms of the parameters V and 0 . This is discussed frcan the 
theoretical point of view in preceding sections. In subsequent 
sections, actual methods of analysis will be described. 

All expansion waves incident upon the wall upstream of the point 
where 0 + v = should be reflected and those incident downstream 
of this point should be neutralized. Thus, this point becomes the 
inflection point of the wall. 

It is interesting to note that, while the initial curve is 
arbitrary, thd corresponding terminal curve is unique once the 
initial curve is established. 

For an infinitely fine mesh of expansion waves, this method is 
exact. Moreover, for a finite mesh size, the finer it is, the more 
accurate is the analysis. 

This method is, perhaps, most useful in designing nonconventional 
nozzles, since for conventional types, the Foelsch method (to be 
described later) is more convenient. 


Puckett’s Method 

Puckett, in reference 7* introduced a variation of Busemann's 
method for designing nozzles. Its advantages will be discussed 
subsequently. The method consists briefly of starting at the middle 
of the nozzle and working toward both ends. 
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The flow through tha nozzle at the maximum expansion section 
(inflection point) is assumed to have a uniform speed and uniformly 
varying direction of flow. Such conditions are illustrated in figure 11. 
The stipulation of these "boundary conditions has "been found from 
experience to "be reasonable. With these boundary values, the term inal 
section of the nozzle can be determined by the same method as for the 
original type Busemarm nozzles. By working backward in a like manner, 
an initial section can also be constructed. Moreover, if 9 1 is less 

than 0 lmax , then one or more of the expansion waves must have been 

reflected. Since there is a choice as to which wave is reflected, 
there is more than one initial curve that will provide the specified 
flow at the maximum expansion section. In fact, if the mesh size 1b 
allowed to became infinitely fine, then it follows that there are an 
infinite number of initial curves that correspond to this terminal 
curve. This same agreement obviously holds for initial curves 
corresponding to other terminal curves. 

While, however, there are an infinite number of suitable initial 
curves for each terminal curve, this does not infer that any contour 
satisfying the area— ratio requirement is suitable . On the other hand, 
the error introduced by using an arbitrary curve can be ignored for 
most practical, purposes, provided that a certain amount of care is 
taken. In a later section a simple method for the design of such 
initial sections will be discussed. 

There are several advantages to Puckett’s method. First, if 
the simplified method for designing the initial section is used, the 
time or work involved in designing a nozzle is approximately halved. 

The second advantage becomes apparent during the actual calcu- 
lation of nozzles. In the original Busemann method, expansion waves 
are originated at certain points along a smooth initial curve; that 
Is, the spacing of the expansion waves is orderly, although It need 
not be uniform. When a finite mesh size Is used, sometimes expansion 
waves are reflected from the wall at such points as to destroy the 
orderliness of the spacing of the ensuing expansion wave pattern. 

The resulting terminal curve thus acquires slight irregularities. 

These irregularities disappear as the mesh size becomes Infinitely 
fine and, in practice, one usually draws a faired curve through them. 

The Puckett method does not avoid this difficulty, but neglects it 
by assuming that the terminal curve is not affected by the wave pattern 
of the initial curve. 


Foelsch’s Method 


Foelsch’s method (reference 3 ) is similar to Puckett’s insofar 
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as one starts at the inflection point of the nozzle and proceeds in 
"both directions. It differs slightly in “boundary conditions “but its 
main difference and direct advantage is that it is analytic. Only 
the portion of Foelsch's theory which deals with the expansion section 
will he discussed here. 

The assumptions of this method., or rather its “boundary conditions, 
may he variously stated (fig. 12): (l) Along the Mach line emanating 

from the inflection point, the velocity vectors are co—original, (2) 
the Mach number is constant along the arc of the circle which passes 
through the inflection point of the wall perpendicularly (and obviously 
its center is the origin of the velocity vectors), (3) in the region 
between this arc and the Mach line from the inflection point, the Mach 
number is a function solely of the radius from the vector origin. 

Using the following notation (fig. 12) 

r distance from vector origin to arbitrary point on inflection 

point Mach line 

r Q hypothetical r for M = 1 

l length of Mach line between inflection-point Mach line and 

final curve 

x coordinates measured from sonic section 

y coordinates measured from center line 

Xi,yi coordinates of inflection point 

Xs,y2 running coordinates of inflection-point Mach line 
y Q sem i height of sonic section of nozzle 
H height of test section (2y^) 

it can be shown that 


r o = (®i in radians) 


r = v r 


(a/a*) 


r i = r Q (Ai/A*) = 


7 1 

sin ©i 


(7) 

(8) 

(9a) 
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or 


2i - (Sl radlans) 

y c 0! \A*/ 

2 = Mr (v — V-jJ (v in radians) 

X2 — xi = — r x cos 0 X + r cos (v^ — v) 

72 = r sin (v t - v) 

and the coordinates of the terminal curve are 

x — Xi=X£ — Xi + 2 cos (vt - V + a) 

y = J2 + 2 sin (v t - V + a) 
the length of the terminal section (in test— section heights) 1 b 

3 2 ~ 1 + 2 ll COS 01 ) ^ 01 ln radlans ) ( x 3 ) 

By varying the Mach number M along the terminal, curve from M x 
to Mfc and determining the corresponding values of a, V, r, and 2, 
the coordinates of the terminal curve can he found and are determined 
as a function of conditions in the test section and at the inflection 
point. Table I is included to facilitate these calculations. The 
initial curve, as for the Puckett method, nay he treated separately. 

It can he seen that the methods of Puckett and Foelsch are quite 
similar with regard to houndary conditions, the former having a constant 
Mach number along a straight line and the latter along a circular arc. 
Both assumptions are eq ually plausible. The difference between these 
assumptions manifests itself in a slight and inconsequential lengthening 
of the Foelsch nozzle relative to the corresponding Puckett nozzle. 

The analytic nature of the Foelsch method allows nozzles to he 
determined to any desired degree of accuracy and without any such 
apparatus (to he described later) as that required for the graphical 
methods. 


(9b) 

( 10 ) 

(11a) 

(lib) 

(12a) 

(12b) 
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It is interestin g to observe that this method is one of the few 
exact analytic solutions of the general nonlinear potential equation 
for compressible supersonic flow. 

The Initial Curve 

The initial curve, either exact or approximate, must satisfy 
certain geometric boundary conditions: It must satisfy the area- 

ratio req uir ements - It must have zero slope at the sonic section 
H.nri the s ame ordinate, slope, and curvature (zero) at the inflection 
point as the terminal curve. The ordinate, slope, and curvature 
should vary monotonically between the sonic section and the inflection 
point. A simple function satisfying these limitations is 

r-To + t <*> 

from which it follows that 

x i = | (7i “ 7 0 ) cot 0 i (15) 


Experience has shown that this approximate initial— curve function 
can be used for both the Puckett and Eoelsch methods without any 
serious error. For the original type Busemann method, this curve can 
be simulated by appropriate straight— line segments. In thiB case the 
curve becomes exact. 


Analysis of Hozzles 

The analysis of given nozzles to determine the velocity distri- 
bution in the test section and ascertain the existence or nonexistence 
of shock waves is a process very similar to the design of nozzles. In 
fact, the procedure for the initial portion of the nozzle is identical. 

In the terminal section, instead of neutralizing the expansion 
waves, they are all reflected and compression waves started at appro- 
priate places. For small angular deflections, compression waves may 
be considered simply as negative expansion waves. In practice, wbere 
an expansion wave is incident upon a wall near the position where a 
compression wave (of the same numerical strength, but opposite sign) 
originates, they may be considered to neutralize each other. 

Thus from the coordinates V and 0, the velocity distribution 
in the test section can be found. The location of possible shock 
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waves is indicated “by a region of converging Mach lines (or compression 
waves) — the greater the concentration of converging waves , the stronger 
the shock. A weak concentration of slowly converging waves may he too 
weak to show up in a schlieren photograph or to have any noticeable 
effect; hence, the term "possible" shock waves ms used. A nozzle 
exhibiting a region of converging waves is shown schematically in 
figure 13. 


EFFECT OF VARIATION OF PARAMETERS 


The major parameter involved in the design of nozzles is 9 ± 
or, perhaps, rather 9 x/— Vj.. The length of the nozzle is intimately 
associated with this parameter. 


As previously stated, the maximum value that 9 ± can have for 
shock— free, flow at a given Mach number is A nozzle so designed 

will be the shortest possible for that Mach number and must have a 


sharp throat such as the one shown in figure 10 with v 


= Iv, 


The 


other extreme in designing nozzles is setting 
require that the nozzle have infinite length. 


0 i = 0 . This would 


There are, of course, certain obvious disadvantages in designing 
a nozzle too long or too short. An excessively long nozzle incurs 
adverse boundary— layer effects of two kinds: First, the longer the 

nozzle, the thicker the boundary layer, other conditions being the 
same. Since boundary— layer thickness is, at the present time, not very 
amenable to precise calculation, a given percent error in boundary- 
layer calculation is more serious when the boundary layer is thick. The 
result is that flow in the test section is less likely to be uniform, 
parallel, and shock free. Second, a thicker boundary layer represents 
an unnecessary mate of energy. 


An excessively short nozzle, on the other hand, is liable to 
other troubles. A minimum— length nozzle has for a given Mach number, 
the maximum number of expansion waves (considering each to be of 
finite strength) concentrated into the minimum space. A longer nozzle 
achieves the same Mach number by reflecting some of the waves. Thus, 
it has fewer of them and these are extended over a wider range. This 
is to say that the expansion waves are more concentrated in shorter 
nozzles. It is then apparent that they are more sensitive to proper 
design than longer ones. Designing nozzles to be somewhat longer than 
the minimum incorporates what might be termed a safety factor. In 
addition, there is less likelihood for such a nozzle to have oscil- 
latory flow. 
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The tendency at some German laboratories was to design nozzles 
with lengths equal to or slightly greater than the minimum. While 
most of these nozzles were claimed to he satisfactory, subsequent 
experience has shown that small gradients previously believed negli- 
gible have been found to exert strong influences on test results. 


Puckett, in reference suggests using Si equal to from one- 
half to two-thirds of Si (Si = — v^) . It is believed that at 


max 2 

low Mach numbers such nozzles will be unnecessarily long. 


At the present time there are insufficient experimental data to 
say exactly how a nozzle should be designed. However, experience up to 
the present time indicates that a value of 


s/s 

Si ~ (Jn) (jr) ( for a±r ) ( i6 ) 

will provide a good working hypothesis for Mach numbers up to about 
five. 


The preceding equation is restricted to air only because of the 
limitations of past experience. The general considerations discussed 
herein, however, apply to helium or any other compressible fluid. 


CONSTRUCTION CG? FLOW FIELD - SUPERSONIC PROTRACTOR 

The determination of the flow field in a nozzle has been discussed 
previously from the theoretical point of view. It remains now to show 
how to construct or determine the orientation of each of the Mach line 
(or expansion wave) segments which make up the net that determines the 
flow field. (See figs. 8 and 9») 

Various methods have been proposed to do this. Analytic methods, 
such as the one described In reference 8, have been devised but are 
extremely tedious. Graphical methods have been found sufficiently 
accurate for most design or anlysis purposes. On the other hand, the 
analytic nature of the Foelsch method allows ordinates to be determined 
simply and precisely. The main use of the Busemann theory is, at the 
present time, usually restricted to the design of nonconventional 
nozzle shapes and the analysis of any given shape. 

A graphical method based on the use of characteristic theory and 
the hodograph plane is described in reference 2. However, this method 
has been superseded by the so-called "supersonic protractor" (reference 9 ) , 
a modification of which is described herein. 
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It is assumed that conditions along an expansion wave are the 
average of those in the regions it separates. Each segment is thus 
characterized by the pair of coordinates V and 9. For each value 
of v and 6 , there are two possible orientations of an expansion 
wave. These correspond to the two Mach lines produced by a point 
disturbance. The angle made by an expansion wave with the datum .line 
is 9 + a for the wave directed upward in the stream direction and 
9 — a for the one directed similarly downward. These two cases are 
shown in figure lk. 

The supersonio protractor has two essential parts which may be 
described independently. The first., shown in figures 15(a) and 15(b) 
consists of a semicircular transparent disk, pivoted at its center, 
and with a straight edge attached. It is graduated along its 
circumference such that when the desired V is set over the datum 
line, for example, v = 3O 0 , a. Is represented as shown. That this 
is possible follows from equation (2) : 

v o°-a) 

The second piece, shown in figure 1 6, consists simply of a 
circular disk graduated along its circumference in degrees. This 
scale represents the stream direction 9. 

If the former part of the protractor, that providing a, is 
rotated through an angle equal, to the stream direction 9, the required 
orientation of the Mach line (or expansion wave) is thus determined. 

This is accomplished with the protractor by superimposing the former 
upon the latter concentrically and rotating the former until the 
desired v is set over the desired 9. This is shown schematically 
in figures 17(a) and 17(b), and the similarity of these with figure 14 
should be noted. Thus, while the end points of certain expansion-wave 
segments may be dependent upon the previous one, each in its turn can 
be orientated simply by means of this protractor, knowing, of course, 
v and 9 . 

Table H, containing values of a corresponding to even values 
of v , is included for calibrating the supersonic protractor. It 
should be noted that if the amount of work involved does not justify 
the construction of this protractor, a drafting machine may be substi- 
tuted. In this case, 9 a can be set with the aid of table H. 



B ODUDAE Y-LAXER CCKKECTIOn 

The problem of correcting nozzle contours for the constricting 
effect of the boundary layer has not been given adequate consideration 
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In the past. However, reference 10, which contains a method for 
est ima t ing boundary— layer thickness, has been found useful. A small 
amount of experimental boundary— layer data is also included in 
reference 7. 

The manner of applying the correction itself is quite simple. 
One merely increases the nozzle ord i nates by an amount equal, to the 
displacement thickness of the boundary layer. It is sometimes 
desirable to keep two walls parallel. In this case, the other two 
walls are corrected to allow for the boundary layer on all four. 


CONCLUDING REMARKS 

Using the methods discussed in this report, it is possible to 
design satisfactory nozzles either graphically or analytically. 
While the analytic method is to be preferred in design, the graphic 
method can be extended to include the analysis of given nozzle 
shapes to determine flow characteristics. A supersonic protractor 
which permits rapid graphical analysis and design is described. No 
correction for boundary layer has been included. 

Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 

Moffett Field, Calif. 


REFERENCES 

1. Busemann, A. : Gasdynamik. Handbuch der Experimental-physik, 

Akademische Verlagsgesellschaft, Leipzig. Bd. 4, Teil 1, 1931* 

S. 34l->t60. 

2. Preiswerk, Erast: Application of the Methods of Gas Dynamics to 

Water Flows with Free Surface. Part I. Flows With No Energy 
Dissipation. NACA TM No. 93^# 19^0° 

3. Foelsch, Kuno: A New Method of Designing Two-Dimensional Laval 

Nozzles for a Parallel and Uniform Jet. North American 
Aviation Rep. No. NA— 46-235, Mar. 1946. 

If. Shapiro, Ascher H.: Nozzles for Supersonic Flow Without Shock 

Fronts. Journal of Applied Mechanics, vol. 11, no. 2, June 1944, 
pp„ A-93 — A— 100. 



1 6 


MCA TTT No. 1651 


5. Liepmann, Hans W., and Puckett , Allen E.: Introduction to 

Aerodynamics of a Compressible Eluid. John Wiley and Sons, 

Inc., 1947 . 

6. Taylor, G. I., and Jfeccoll, J. W. : The Mechanics of Compressible 

Eluids. Two-Dimensional Flow at Supersonic Speeds. Aerodynamic 
Theory, vol. HI, div. H, ch. IV, W. F. Durand, ed., 

Julius Springer (Berlin), 1935, PP. 234—249. 

7. Puckett, Allen E.: Supersonic Nozzle Design. Journal of Applied 

Mechanics, vol. 13, no. 4, Dec. 1946, pp. A-265 — A-27O. 

8. Shapiro, Ascher H., and Edelman, Gilbert M. : Method of Character- 

istics for Two-Dimensional Supersonic Flow — Graphical and 
Numerical Procedures. Journal, of Applied Mechanics, vol. 14, 
no. 2, June 1947, pp. A-154 — A— 162. 

9. Temple, G. ; The Method of Characteristics in Supersonic Flow. 

British Rep. No. S.M.E. 3^82, R.A.E., Jan. 1944-. 

10. Puckett, Allen E.; Final Report on the Model Supersonic Wind- 

Tunnel Project. NDRC, Div. 2, Armor and Ordnance Rep. A-2 69. 
(Office cf Scientific Research and Development) Rep. 3569, 

Apr. 30, 1944. 



IS 


MCA TIT Ho. 1651 


IT 


TABLE I.- ESSENTIAL PARAHETffiS USED II BQZZLS DESIGN 1 
[ 7 - 1.400 (air) J 


M 

p 

Po 

A* 

A 

a 

(deg) 

V 

(deg) 

M 

P 

Po 

A* 

k ~ 


V 

(deg) 

1.00 

0.5283 

1.0000 

90.00 

0 

1.50 

0.2724 

0.8502 

4l.8i 

11.91 

1.01 

.5221 

. 9999 

81.93 

.04473 

1.52 

.2646 

.84o4 

4l.l4 

12.49 

1.02 

.5160 

.9997 

78.64 

.1257 

1.54 

.2570 

.8304 

4o.49 

13.09 

1.03 

.5099 

.9993 

76.14 

.2294 

1.56 

.2496 

.8203 

39.87 

13.68 

1.04 

.5039 

.9987 

74.06 

.3510 

1.53 

.2423 

.8101 

39-27 

14.27 

1.05 

.4979 

.9980 

72.25 

.4874 

1.60 

.2353 

.7998 

38.68 

14.86 

1.06 

.4919 

.9971 

70.63 

.6367 

1.62 

.2284 

.7895 

38.12 

15.45 

1.07 

.4860 

.9961 

69.16 

.7973 

1.64 

.2217 

•7791 

37.57 

16. 04 

1.08 

.4800 

.9949 

67.81 

.9680 

1.66 

.2151 

.7686 

37.04 

16.63 

1.09 

.4742 

•9936 

66.55 

1.148 

1.68 

.2088 

.7581 

36.53 

17.22 

1.10 

.4684 

.9921 

65.38 

1.336 

1.70 

.2026 

.7476 

36.03 

17.81 

1.11 

.4626 

.9905 

64.28 

1.532 

1.72 

.1966 

.7371 

35.55 

l8.4o 

1.12 

.4568 

.9888 

63.23 

1.735 

1.74 

.1907 

.7265 

35.08 

18.98 

1.13 

.4511 

.9870 

62.25 

1.944 

1.76 

.1850 

.7160 

34.62 

19.56 

1.14 

.4455 

.9850 

61.31 

2.160 

1.78 

.1794 

.7054 

34.18 

20.15 

1.15 

.4398 

.9828 

60. 4i 

2.381 

1.80 

.1740 

.6949 

33.75 

20.73 

1.16 

.4343 

.9806 

59.55 

2.607 

1.82 

.1683 

.6845 

33.33 

21.30 

1.17 

.4287 

.9782 

53.73 

2.839 

1.84 

.1637 

.6740 

32.92 

21.88 

1.18 

.4232 

•9758 

57.94 

3.074 

1.86 

.1587 

.6636 

32.52 

22.45 

1.19 

.4178 

•9732 

57.18 

3.314 

1.88 

.1539 

.6533 

32.13 

23.02 

1.20 

.4124 

.9705 

56.44 

3.558 

1.90 

.1492 

.6430 

31.76 

23.59 

1.21 

.4070 

.9676 

55.74 

3.806 

1.92 

.1447 

.6328 

31.39 

24.15 

1.22 

.4017 

.9647 

55.05 

4.057 

1.94 

.1403 

.6226 

31.03 

24.71 

1.23 

.'3964 

.9617 

54.39 

4.312 

1.96 

.1360 

.6125 

30.68 

25.27 

1.24 

.3912 

.9586 

53.75 

4.569 

1.98 

.1318 

.6025 

30.33 

25.83 

1.25 

.3861 

.9553 

53.13 

4.830 

2.00 

.1278 

.5926 

30.00 

26.38 

1.26 

.3809 

.9520 

52.53 

5.093 

2.02 

.1239 

.5328 

29.67 

26.93 

1.27 

.3759 

.9486 

51.94 

5.359 

2.04 

.1201 

.5730 

29.35 

27.48 

1.28 

.3708 

.9451 

51.38 

5.627 

2.06 

.1164 

.5634 

29-04 

28.02 

1.29 

.3658 

.9415 

50.82 

5.898 

2.08 

.1128 

.5538 

28.74 

28.56 

1.30 

.3609 

.9378 

50.28 

6.170 

2.10 

.1094 

.5444 

28.44 

29.10 

1.31 

.3560 

.9341 

49.76 

6.445 

2.12 

.1060 

.5350 

28.14 

29.63 

1.32 

.3512 

.9302 

49.25 

6.721 

2.14 

.1027 

.5258 

27.86 

30.16 

1.33 

.3464 

.9263 

48.75 

7.000 

2.16 

.09956 

.5167 

27.58 

30.69 

1.34 

.3417 

.9223 

48.27 

7.279 

2.18 

.09650 

.5077 

27.30 

31.21 

1.35 

.3370 

.9182 

47.79 

7.561 

2.20 

.09352 

.4988 

27-04 

31.73 

1.36 

.3323 

.9141 

47.33 

7.844 

2.22 

.09064 

.4900 

26. T7 

32.25 

1.37 

.3277 

.9099 

46.83 

8.128 

2.24 

.08785 

.4813 

26.51 

32.76 

1.38 

.3232 

.9056 

46.44 

8.413 

2.26 

.08514 

.4727 

26.26 

33.27 

1.39 

.3187 

.9013 

46.01 

8.699 

2.28 

.08252 

.4643 

26.01 

33.78 

l.4o 

.3142 

.8969 

45.58 

8.987 

2.30 

.07997 

.4560 

25.77 

34.28 

1.42 

.3055 

.8880 

44.77 

9.565 

2.32 

.07751 

.4478 

25.53 

34.78 

1.44 

.2969 

.8788 

43.98 

10.15 

2.34 

.07512 

.4397 

25.30 

35.28 

1.46 

.2886 

.8695 

43.23 

10.73 

2.36 

.07281 

.4317 

25.07 

35.77 

1.48 

.2804 

• 8599 

42.51 

11.32 

2.38 

.07057 

.4239 

24.85 

36.26 


Adapted from Notes and Tables for Use in the Analysis of Supersonic Plov 
by the staff of the Ames 1— by 3-foot supersonic vind— tunnel section, 
MCA TN Ho. 1428, 1947 . 
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TABLE I.- CONTINUED. ESSENTIAL PARAMETERS USED IN NOZZLE DESIGN 1 


M 

L_ 

Po 

£1 

A 

a 

(deg) 

V 

(deg) 

M 

E-xlO 3 

Po 

A 

(deg) 

V 

(deg) 

2. 40 

0.06840 

0.4161 

24.62 

36.75 

4.00 

6.586 

0.09329 

14.48 

65.78 

2.42 

.06630 

.4085 

24. 4l 

37.23 

4.10 

5.769 

.08536 

14.12 

67.08 

2.44 

.06426 

.4010 

24.19 

37-71 

4.20 

5.062 

.07818 

13.77 

68.33 

2.46 

.06229 

.3937 

23.99 

38.18 

4.30 

4.449 

.07166 

13.45 

69.54 

2.48 

.06033 

.3864 

23.78 

38.66 

4.40 

3.918 

.06575 

13.14 

70.71 

2.50 

•05853 

.3793 

23.58 

39-12 

4.50 

3.455 

.06038 

12.84 

71.83 

2.52 

.05674 

.3722 

23.38 

39.59 

4.60 

3.053 

.05550 

12.56 

72.92 

2.54 

.05500 

• 3653 

23.18 

40.05 

4.70 

2.701 

.05107 

12.28 

73.97 

2.56 

.05332 

• 3585 

22.99 

40.51 

4.80 

2.394 

.04703 

12.02 

74.99 

2.58 

.05169 

• 3519 

22.81 

40.96 

4.90 

2.126 

•04335 

11.78 

75.97 

2.60 

.05012 

.3453 

22.62 

4l.4l 

5.00 

1.890 

.04000 

11.54 

76.92 

2.62 

.04859 

.3389 

22.44 

41.86 

5.1 

1.683 

.03694 

11.31 

77.84 

2.64 

.04711 

.3325 

22.26 

42.31 

5.2 

1.501 

.03415 

11.09 

78.73 

2.66 

.04568 

.3263 

22.08 

42.75 

5-3 

1.341 

.03160 

10,88 

79-60 

2.68 

.04429 

.3202 

21.91 

43.19 

5.4 

1.200 

.02926 

10.67 

80.43 

2.70 

.04295 

.3142 

21.74 

43.62 

5.5 

1.075 

.02712 

10.43 

81.24 

2.72 

.04165 

.3083 

21.57 

44.05 

5.6 

.9643 

.02516 

10.29 

82.03 

2.74 

.04039 

.3025 

21.41 

44.48 

5.7 

.8664 ■ 

.02337 

10.10 

82.80 

2.76 

.03917 

.2968 

21.24 

44.91 

5.8 

.7794 

.02172 

9.928 

83.54 

2.78 

•03799 

.2912 

21.08 

45.33 

5.9 

.7021 

.02020 

9.758 

84.26 

2.80 

.03685 

.2857 

20.92 

45.75 

6.0 

.6334 

.01880 

9.594 

84.96 

2.82 

.03574 

.2803 

20.77 

46.16 

6.1 

• 5721 

.01752 

9.435 

85.63 

2.84 

.03467 

.2750 

20.62 

46.57 

6.2 

.5174 

.01634 

9.282 

86.29 

2.86 

-03363 

.2698 

20.47 

46.98 

6.3 

,4684 

.01525 

9.133 

36.94 

2.88 

.03263 

.2648 

20.32 

47.39 

6.4 

.4247 

.01424 

8.989 

87.56 

2.90 

.03165 

.2598 

20.17 

47.79 

6.5 

.3855 

.01331 

8.850 

88.17 

2-92 

.03071 

.2549 

20.03 

48.19 

6.6 

.3503 

.01245 

8.715 

88.76 

2.94 

.02980 

.2500 

19.89 

43.59 

6.7 

.3187 

.01165 

8.584 

•39.33 

2.96 

.02891 

.2453 

19.75 

48.98 

6.8 

.2902 

.01092 

8.457 

89.89 

2.98 

.02805 

.2407 

19. 6i 

49.37 

6.9 

.2646 

.01024 

8.333 

90.44 

3.00 

.02722 

.2362 

19.47 

49.76 

7.0 

.2416 

.009602 

8.213 

90.97 

3.05 

.02526 

.2252 

19.14 

50.71 

7.1 

.2207 

.009015 

8.097 

91.49 

3-10 

.02345 

.2147 

18.82 

51.65 

7-2 

.2019 

.008469 

7.984 

92.00 

3.15 

.02177 

.2048 

18.51 

52.57 

7.3 

.1848 

.007961 

7.873 

92.49 

3.20 

.02023 

.1953 

18.21 

53.47 

7.4 

.1694 

.007490 

7.766 

92.97 

3.25 

.01880 

.1863 

17.92 

54.35 

7-5 

.1554 

.007050 

7.662 

93.44 

3.30 

.01748 

.1777 

17.64 

55.22 

7.6 

.1427 

.006641 

7.561 

93.90 

3-35 

.01625 

.1695 

17.37 

56.07 

7.7 

.1312 

.006259 

7.462 

94.34 

3.40 

.01513 

.1617 

17.10 

56.91 

7.8 

.1207 

.005903 

7.366 

94.76 

3-45 

.Ol4o8 

.1543 

16.85 

57-73 

7-9 

.1111 

.005571 

7.272 

95.21 

3-50 

.01311 

.1473 

16.60 

58.53 

8.0 

.1024 

.005260 

7.181 

95.62 

3-60 

.01138 

.1342 

16.13 

60.09 

8.1 

.09448 

.004970 

7.092 

96.03 

3.70 

.009903 

.1224 

15.68 

61.60 

8.2 1 

.08723 

.004698 

7.005 

96.43 

3.80 

.008629 

.1117 

15.26 

63.04 

8.3 

.08060 

.004444 

6.920 

96.82 

3.90 

.007532 

.1021 

14.86 

64.44 

8.4 

.07454 

.004206 

6.837 

97.20 


1 Adapted from notes and Tables for Use in the Analysis of Supersonic Flo* by 
the staff of the Ames 1— by 3— foot supersonic wind-tunnel section. NACA 
TIT No. 1428, 1947. 
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TABLE I.- CONCLUDED. 

ESSENTIAL PARAMETERS USED IN NOZZLE DESIGN 1 


M 

— XlO 3 

Po 

— XlO 3 
A 

a 

(deg) 

V 

(leg) 

8.5 

0.06896 

3.981 

6.756 

97-58 

8.6 

.06390 

3.773 

6.677 

97.94 

8.7 

.05923 

3.577 

6.600 

98.29 

8.8 

. 05494 - 

3.392 

6.525 

98.64 

8.9 

.05101 

3.219 

6.451 

98.98 

9..0 

.04739 

3.056 

6.379 

99.32 

9-1 

.04405 

2.903 

6.309 

99.65 

9.2 

.04099 

2.759 

6.240 

99.97 

9.3 

.03816 

2.623 

6.173 

100.28 

9.k 

.03555 

2.495 

6.107 

100.59 

9-5 

.03314- 

2.374- 

6.042 

IOO .89 

9-6 

.03092 

2.261 

5-979 

101.19 

9-7 

.02886 

2.153 

5.917 

101.48 

9.8 

.02696 

2.052 

5.857 

101.76 

9-9 

.02520 

1.956 

5-797 

102.04- 

10.0 

.02356 

1.866 

5-739 

102.32 


'’'Adapted from Notes and Tables for Use in 
the Analysis of Supersonic Flow by the 
staff of the Ames 1— by 3— foot supersonic 
wind-tunnel section. MCA TN No. 1428^ 
19 ^ 7 . 
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TABLE II.- PARAMETERS USED IE CALIBRATIHG SUPERSONIC PROTRACTOR 1 

[7 = 1.400 (air)] 


V 

(deg) 

M 

a 

(deg) 

V 

(deg) 

M 

a 

(deg) 

0 

1.0000 

90.00 

44 

2.7179 

21.59 

1 

1.0808 

67.70 

45 

2.7643 

21.21 

2 

1.1328 

61.96 

46 

2.8120 

20.33 

3 

1.1770 

58.17 

47 

2.8610 

20.46 

4 

1.2170 

55.29 

48 

2.9105 

20.09 

5 

1.2554 

52.77 

49 

2.9616 

19.73 

6 

1.2935 

50.63 

50 

3.0131 

19.38 

7 

1.3300 

48.75 

51 

3.0660 

19.06 

8 

1.3649 

47.11 

52 

3.U93 

18„70 

9 

1.4005 

45.57 

53 

3.1737 

18.38 

10 

1.4350 

44.18 

54 

3.2293 

18. 04 

ll • 

1.4688 

42.92 

55 

3.2865 

17.72 

12 

1.5028 

41.72 

56 

3.3^51 

17.40 

13 

1.5365 

40.60 

57 

3.4055 

17.08 

14 

1.5710 

39.53 

58 

3.4675 

l5„76 

15 

1.6o45 

38.54 

59 

3.5295 

160 46 

16 

1.6380 

37.63 

60 

3.5937 

16.16 

17 

1.6723 

36.73 

62 

3.7283 

15.56 

18 

1.7061 

35.88 

64 

3.8690 

14.98 

19 

1.7401 

35.08 

66 

'+.0164 

14.42 

20 

1.7743 

34.31 

68 

4.1733 

13.86 

21 

1.8090 

33.54 

70 

4.3385 

13.33 

22 

1.8445 

32.83 

72 

4.5158 

12.79 

23 

I08795 

32.15 

74 

4.7031 

12.28 

24 

1.9150 

31.49 

76 

4.9032 

11.76 

25 

1.9502 

30.85 

78 

5.U9 

11.27 

26 

1.9861 

30.23 

80 

5.349 

10.78 

27 

2.0222 

29.64 

82 

5.595 

10.29 

28 

2.0585 

29.06 

84 

5.867 

9.81 

29 

2.0957 

28.49 

86 

6.155 

9.35 

30 

2.1336 

27.97 

88 

6.472 

8.88 

31 

2.1723 

27.41 

90 

6.820 

8.43 

32 

2.2105 

26.90 

92 

7.202 

7.98 

33 

■ 2.2492 

26.40 

94 

7.623 

7-54 

3^ 

2.2885 

25.91 

96 

8.093 

7.10 

35 

2.3288 

25.43 

98 

8.622 

6.67 

36 

2.3688 

24.99 

100 

9.210 

6.23 

37 

2.4108 

24.53 

102 

9.887 

5.80 

38 

2.4525 

24.07 

104 

10.658 

5.38 

39 

2.4942 

23.64 

108 

12.58 

4.56 

40 

2.5372 

23.22 

112 

15.37 

3.73 

4l 

2.5810 

22.80 

116 

19.70 

2.91 

42 

2.6254 

22.38 

120 

27.29 

2.10 

4-3 

2.6716 

21.98 

124 

44.08 

1.30 


^Adapted from reference 7„ 
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Hgure. I- Section of supersonic wind funnel showing relative 
position of nozzle. 
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Figure 2.- Supersonic flow about c* comer: 



Figure 3- Supersonic flow about a comer 
from an initial Mach number of unify . 
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Figure, 4 - Approximate representation 
of flow about a small corner. 



Figure, 6.- Intersection of two expansion 
Waves. 
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Figure 3. - Flow fiofd in a /7oz2/e with no Wares reffeohed. 
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Fiaure 13.- Nozzle <3xh}b/Hnq a reqion of conve-nq/nq 

com/On&s&ion) waves and uncan c ei/ed expansion waves. 
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Figure 14-,- Mach lines from a. point 
disturbance. 
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(cl) Mach line directed upward 
in stream direction 

Stream 



(b) Mach tine directed downward 
in stream direction 


HguLre U5- Upper half of supersonic 
protractor . 
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strea m cUrec f/'o n 



(b) Mach line directed downward 
in stream direct/ on 

Figure f7- Assembly of supersonic protractor. 


HACA - Langley Field, Vj. 



